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Abstract 

This work investigates the physics of elementary excitations for the 
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so-called relativistic quantum scalar plasma system, also known as the 
Higgs- Yukawa system. Following the Nemes-Piza-Kerman-Lin many- 
body procedure, the Random-Phase Approximation (RPA) equations 
were obtained for this model by linearizing the Time-Dependent Hartree- 
Fock-Bogoliubov equations of motion around equilibrium. The result- 
ing equations have a closed solution, from which the spectrum of ex- 
citation modes are studied. We show that the RPA oscillatory modes 
give the one-boson and two-fermion states of the theory. The results 
indicate the existence of bound states in certain regions in the phase 
diagram. Applying these results to recent LHC observations concern- 
ing the mass of the Higgs boson, we determine limits for the intensity 
of the coupling constant g of the Higgs- Yukawa model, in the RPA 
mean-field approximation, for three decay channels of the Higgs bo- 
son. Finally, we verify that, within our approximations, only Higgs 
bosons with masses larger than 190 GeV/c 2 can decay into top quarks. 
PACS Numbers: ll.10.Lm, ll.10.Gh, 03.65.Nk, 03.65. Ge, 21.60.Jz 

Keywords: Higgs- Yukawa system. Nemes-Piza-Kerman-Lin proce- 
dure. Random-Phase Approximation. Bound states. 

1 Introduction 

Recent years have witnessed substantial progress towards understanding the 
non-equilibrium time evolution of quantum fields. Results have been ob- 
tained that proved important in several applications. Examples are found 
in cosmology, such as the description of quantum-field expectation values in 
the early universe and subsequent hot stage (big bang) [1, 2], in high en- 
ergy particle physics, such as the description of the dynamics in heavy-ion 
collision experiments, which seeks to establish experimental signatures for 
the non-equilibrium evolution of the quark-gluon system and chiral phase 
transition [3, 4, 5], and in complex many-body quantum systems, such as the 
description of the dynamics of the Bose Einstein condensates [6, 7], among 
other applications [8]. 

In a previous publication we have developed a framework to investigate 
the initial-value problem in the context of interacting fermion- scalar field 
theories [9]. This framework had earlier been developed in the context of the 
non-relativistic nuclear many-body dynamics by Nemes and Toledo Piza [10]. 



The method led to a set of self-consistent time-dependent equations for the 
expectation values of linear and bilinear forms of field operators. These dy- 
namical equations acquire a kinetic type structure in which the lowest-order 
approximation corresponds to the usual Gaussian mean-field approximation. 
As an application, a zero-order calculation was implemented within the sim- 
plest context of a relativistic quantum scalar plasma system, nowadays also 
called the Higgs- Yukawa system. The usual renormalization prescription was 
shown to be also applicable to this non-perturbative calculation. In partic- 
ular, a finite expression for the energy density was obtained, the numerical 
results suggesting that the system always has a single stable minimum. 

Here we report an application of the renormalized nonlinear dynamical 
equations obtained in our previous article [9] and follow Kerman and Lin 
[11, 12] to investigate the near-equilibrium dynamics around the stationary 
solution of a Higgs- Yukawa system. We will show that one-boson and two- 
fermion physics can be studied in the linear approximation of the mean-field 
equations. We will solve those equations in closed form and find the scattering 
amplitude, as well as the conditions allowing a two-fermion bound state. 

The motivation for this work is the recent observation of a possible Higgs 
boson around 125 Gev/c 2 [13]-[15]. Since the early 1990's, the Higgs- Yukawa 
model has been used to better understand the fermion mass generation via 
the Higgs mechanism. Recently, the Higgs- Yukawa model has been used to 
impose limits on the Higgs mass and on the intensity of the Higgs- Yukawa 
coupling. It has also been used to study the consequences of heavy extra- 
generation fermions (mass > 600/c 2 GeV). An important consequence of a 
fourth fermion generation is the possibility of formation of bound states that 
can replace the role of the Higgs boson. Recently, these issues have been 
intensively studied on the lattice [16] or perturbatively — the 1/iV-expansion 
[17], for example. Our work carries out a non-perturbative calculation of the 
ground state (vacuum) for an interacting Higgs- Yukawa system. 

For clarity and notational purposes, a few key equations from [9] are re- 
peated here. The dynamics of the relativistic quantum scalar plasma model, 
or Higgs- Yukawa model, is governed by the following lagrangian density [18]- 
[30]: 

L = ipli'-f.d — m)ip + gip<f)ip 
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The model Hamiltonian is given by the equality 
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with the shorthand / x = / d 3 x. 

In the Heisenberg picture, 4>(x) and U(x) are scalar spin-0 fields expanded 

as 
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where fcl(t) and b p (t) are boson creation and annihilation operators, po = 
y/p 2 + Q 2 , f2 is the mass parameter of the bosonic fields, p x = p^t — p.x, 
and ip{x) and ^(x) are fermionic spin-1/2 fields, 
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where o^ (£) and o^^i) [a^ 8 (*) an d a ksW] are fermion creation and anni- 
hilation operators associated with the positive- (negative-) energy solutions 
«i(k, s) [«2(k, s)] of Dirac's equation. Likewise, k = y/k 2 + M 2 , M is the 
mass parameter of the fermionic fields and k x = k$t — k.x . 



In Eq. (2), the parameters m and fi are the masses of the fermion and 
of the scalar particles, respectively, and g is the coupling constant. The last 
term on the-right hand side, which encompasses the counterterms necessary 
to remove the later-occuring infinities, is given by the expression 



AttHc 



A Am 6 ^^ 



+ §^ + ** 4 -f<w 



(5) 



where the coefficients A, 5fi 2 , C and D are infinite constants, to be defined 
later. To study the dynamical RPA regime, we have to introduce the wave- 
function renormalization constant Z [23]. 

To deal with condensate and pairing dynamics of the scalar and fermionic 
fields, we first define the unitary Bogoliubov transformation for the bosonic 
sector as follows [31]: 
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where d p is the shift boson operator defined by the expression 



d p (t) = b p {t) - B p (t)8 pfi with B p (t) = (b p (t)) = Tr BF [b p (t)F] . (7) 

Here, J 7 is the unitary many-body density operator [9] describing the system 
state, and the symbol Tr BF denotes a trace over both bosonic and fermionic 
variables. Partial traces over bosonic or fermionic variables will be denoted 
Tr B and Tr F , respectively. For simplicity, we restrict our treatment of the 
fermionic sector to the Nambu transformation [32], parameterized in the 
following form, which incorporates the unitarity constraints: 
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The procedure adopted in [9] to obtain the equations of motion for the 
Nambu-Bogoliubov parameters (pk(t) , 7k (t), ?7p(£), /%>(£), for the quasi- 

particle occupation numbers v\^' s = (a k j a^J), v p = (/3p/3 p ), and for 
the condensates ((f)) and (II) had been developed earlier in the context of 
the non-relativistic nuclear many-body dynamics by Nemes and Toledo Piza 
[10]. That approach follows the line of thought of a time-dependent pro- 
jection technique proposed by Willis and Picard [33] in the context of the 
master equation for coupled systems. The method consists, essentially, of 
writing the correlation information of the full density of the system J 7 in 
terms of a memory kernel acting on the uncorrelated density J-q, with the 
help of a time- dependent projector. At this point, a systematic mean-field 
expansion for two-point correlations can be performed [10]. The lowest order 
corresponds to the results of the usual Gaussian approximation [32, 34]. The 
higher orders describe the dynamical correlation effects between the sub- 
systems and are expressed by means of suitable memorial integrals added 
to the mean-field dynamical equations. The resulting equations acquire the 
structure of kinetic equations, with the memory integrals playing the role of 
collisional dynamics terms. This systematic expansion scheme for memory 
effects, in which the mean energy is conserved to all orders [35, 36, 37, 38], 
was implemented, for example, for the Jaynes-Cummings system [39]. 

In this context, to study the near-equilibrium dynamics around the sta- 
tionary solution of the Higgs- Yukawa system, Ref. [9] focused on a selected 
set of Gaussian observables, which are related to the expectation values of 
linear, 4>(x), H(x), and bilinear forms of field operators, such as <j)(x)(j)(x) , 
I 4>{x)'^{x) , il}(x)ij)(x), etc.. The time evolution of such quantities obeys the 
Heisenberg equation of motion 

i(6)=Tr w [0,H]F , (9) 

where O can be <f>(x), (f)(x)(f)(x), ip(x)ip(x), etc., and T is the state of the 
system, which is assumed to be spatially uniform, in the Heisenberg picture. 
As an approximation, we replace the full density J 7 by a truncated ansatz 
•7"o(£), whose trace is also unitary, and which implements the double mean- 
field approximation [9]. By construction, Tq is written as the most general 
Hermitian Gaussian functional of the field operators consistent with the as- 
sumed uniformity of the system. It will thus be written as the exponential 
of a general quadratic form in the field operators, which can be reduced to 
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diagonal form by a suitable canonical transformation. In this way, J-q repro- 
duces the corresponding J 7 averages for linear or bilinear field operators [32] . 
In particular, we have used a formulation appropriate for the many-body 
problem, where Tq is written in the momentum basis as [32, 34, 39] 
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where a (a^) and /3 (/r ) stand for Nambu-Bogoliubov quasi-particle annihi- 
lation (creation) operators for fermions and bosons, respectively, z/ k s ( z/ p ) 
are the quasi-fermion (quasi-boson) occupation numbers, and A = 1 (A = 2) 
is associated with the positive- (negative-) energy solutions. 

From Eqs. (9)-(ll) we can directly obtain the equation of motion for 
the occupation numbers, for the Gaussian variables, now represented by the 
Nambu-Bogoliubov parameters, and for the condensates. The following ex- 
pressions result: 
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where we have introduced the notation 
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In Eqs. (12)-(19) the quantities M and f2 are the mass parameters of the 
fermionic ip and bosonic fields of the Hamiltonian (2), while m = m — g(<f)) 
stands for the effective mass of fermion particle [9] . 

Another physical quantity of interest is the mean-field energy density of 
the system, 
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The above equations describe the real-time evolution of the Higgs- Yukawa 
system in the double Gaussian mean-field approximation. The results ob- 
tained in (12)-(20), as discussed in [9], are consistent with those in the lit- 
erature, obtained via different approaches, in particular with those obtained 
in Ref. [22] on the basis of a Vlasov-Hartree approximation. 



Reference [9] showed, in detail, that the usual form of renormalization 
[22] is applicable to the non-perturbative procedure described in that pa- 
per. A simple numerical calculation has also shown that the system always 
has a single stable minimum, although, as it has been suggested [9], addi- 
tional investigation is necessary concerning oscillatory modes. The standard 
approach to this question uses the RPA analysis, the resulting eigenvalues 
giving an indication of stability [11, 12, 23]. 

Finally, we note that dynamical correlation corrections can in principle 
be systematically added to the double- Gaussian mean-field calculations with 
the help of a projection technique discussed in [10, 37, 39]. The occupation 
numbers are then no longer constant, a modification that affects the effec- 
tive dynamics of the Gaussian observables. The framework presented in this 
paper also serves as groundwork for finite-density and finite-temperature dis- 
cussions [40]. In particular, a finite matter-density calculation beyond the 
mean-field approximation allows one to study such collisional observables as 
the transport coefficients. The extension of this procedure to nonuniform 
systems is straightforward, albeit long. In this case, the spatial dependence 
of the field are expanded in natural orbitals of the extended one-body density. 
A more general Bogoliubov transformation [31] would relate these orbitals 
to a momentum expansion. 

The results in Eqs. (12)-(18) are nonlinear time-dependent field equations. 
A closed solution is not easily constructed. Here we consider those equations 
in the small oscillation regime and find a closed solution offering insight into 
divers properties of the theory. 

The paper is structured as follows. In Sec. 2, the RPA equations are 
derived for this model by considering near-equilibrium dynamics around the 
stationary solutions obtained in Ref. [9]. Section 3 finds analytical solutions 
for the RPA equations by using a well-know procedure from the scattering 
theory. Section 4 discusses renormalization within the context of scattering 
amplitudes and discusses the existence of bound state solutions. Finally, we 
apply these results to find limits to the intensity of the coupling constant g 
of the Higgs- Yukawa model, in the RPA mean-field approximation, for three 
decay channels of the Higgs boson. Section 5 presents our conclusions. 



2 Near equilibrium dynamics 

The energy density (20) is a function of the Nambu-Bogoliubov parame- 
ters (pk(t) , 7k(0> Vp(t)j k p(0> °f the quasi-particle occupation numbers 

u ks = { a ki a k,s)i v p = Wpl^p )> an d of the condensates (0) and (II). The 
minimum in Eq. (20) corresponds to the ground state of the system. The 
small amplitude motion around the minimum is obtained by linearization of 
the Gaussian motion equations (12)-(18), yielding a set of harmonic oscilla- 
tors [41]. The eigenvalues and the normal modes of these small oscillations 
are the RPA solutions. Physically, the RPA solutions are seen as the energy 
and the wave functions of quantum particles. This section derives the RPA 
equations of the model, whose solutions are discussed in the Sec. 3. 

First, we have to consider the stationary problem [9]. We only have 
to recall Eqs. (12)-(18) to see that 7k = <pk = k p = f) p = (0) = (II) = 
under stationary conditions. Reference [9] discussed the renormalization 
conditions and the solutions for this set of stationary equations in detail. 
In particular, for the renormalization coefficients of "H c , the following self- 
consistency renormalization prescription was chosen [9, 23]: 

D = ±48irg 4 L(m) (21) 

S/j, 2 = T 247rg 4 L(m)G(/j,)TlQng 2 G(0)±247rm 2 g 2 L(m) (22) 

C = ^487rm^ 3 L(m) (23) 

A = ±24nmg 3 L(m)G(fi) ± 16nmgG(m) (24) 

with 



L(m) "ik2k^T^^ ' (25) 

where M — m, without loss of generality [9]. 

The substitution of such counterterms in the stationary equations yields 
the appropriate cancelations, which makes the equations finite, except for 
the combination of the type L(m)[G(fi) — G(Q)]. Since Q is an arbitrary 
expansion mass parameter, one can remove this divergence by setting Q = \x 
[9]. The resulting finite stationary equations for the system can be regrouped 
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as follows: 
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Equations (26)-(31) can be numerically solved for any given \x and g, in 
units of 777, as shown by Ref. [9]. 

To obtain the near-equilibrium dynamics (RPA regime), we examine the 
fluctuations around the stationary solution, namely, 
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where y3 k | oq , 7 k | cq , 77 p | eq , /%,| cq , (</>)| oq and (II)| oq satisfy Eqs. (26)-(31) and 
the deviations 5<pk, 5j^, 5r] p , 8k p , S{4>) and (5(11) are assumed to be small. 
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Next, we expand Eqs. (13)-(18) to first order in the fluctuations. The 
following equations result: 

|k| 
5<Pk = gi&L-j-hk (33) 

3(0>| cq |k|<57 k = -4ko(k 2 + m 2 )5 Vk -2g\k\k 5( ( f ) } (34) 

m = J^z) 6 ^) (35) 

*<n> = -^ + ^ + ^G(^5(0)-g(0)| eq 5(0) 

In the RPA regime, the bosonic variables show no dynamical evolution. 

To eliminate the quantities 5^ and 5(11), can be eliminated by differen- 
tiating Eqs. (33) and (35) with respect to time, so that Eqs. (33)- (36) are 
rewritten as the second-order differential equations 

5<p k = -4k 2 5 Vk -2g\k\5(<t ) ) (37) 

(1 + Z)6(4>) = -(/i 2 + s)5(0) +167rgfh(k')S(p k/ , (38) 

with the notation 

h(k) = B (39) 

Kq 

where 

k = Vk 2 + m 2 (40) 

S = ^ 2 + ^4 7 rG( ;U ) + C(0)| oq + ^((0)| cq ) 2 . (41) 

The small oscillation-dynamics of the Higgs- Yukawa system is therefore 
described by coupled equations of linear oscillators, as usual in the RPA 
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treatment [41]. In particular, when g = these modes are decoupled and 
yield two equations describing simple oscillators. 

To solve the problem (37)- (38) we have to determine the normal modes 
of the small oscillations and their frequencies. Earlier studies have demon- 
strated that these elementary excitations can be interpreted as quantum 
particles. In our case, 5(pk can be seen as two-fermion spinless wave func- 
tion [42], while 5{<j)) provides the one-boson physics of the system [11]. The 
relative momentum of the two-fermion states is |k| [42], while in the scalar 
sector, the particles have no momentum dependence. 

3 RPA equations as a scattering problem 

Section 2 obtained the linear approximation for the Gaussian equations of 
motion, Eqs. (37)- (38). We will now show that these coupled linear oscillator 
equations can be analytically solved, to determine the wave functions and the 
elementary-excitation spectrum of our system. 

We first consider the Fourier transform of the wave functions in the energy 
representation, i. e., the standard relations 



S(p k (t) = du 8(p k (uj) e iuJt 

(42) 
S((f))(t) = ( duj 5{<P)(uj) e wt , 



where 5(pk(u) and S{4>) (u) are now energy-dependent amplitudes. 

We then substitute Eq. (42) into (37)-(38), to obtain the following equa- 
tions: 



(u 2 - All) ^k(w) = 2g\k\5(<P)(uj) (43) 

(w 2 -/i 2 + Zw 2 -S)5(0)(w) = -lQirg [h(k')5ipv(uj) . (44) 

Since the oscillation amplitudes in Eqs. (43-44) play the roles of wave 
functions of quantum particles, it is more convenient to treat this system as 
a coupled-channel scattering problem with appropriate boundary conditions 
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[12]. The following discussion will focus on the scattering process, where the 
source is a two-fermion wave. In this case, from Eq. (44) we have that 



*«<<"> = L-J^-s ) L fc(k '» ^ ■ (45) 



Substitution of Eq. (45) into (43) then yields the result 
/u; 2 -4fc 2 \ / -32irg 2 \ |k| /■ |k'| . , . 

{-^) 5 ^) = ( h;a .^ + z ^_ E J y 1 ^M«) ' (46) 

where the Green's Function includes the effects of coupling Sip^ to 5 ((f)). 
The potential is separable [43], in the sense that 

(k|T/|k'> = „(k) v(k>) = B H . (47) 

/CO /Cq 

In the general solution of Eq. (46), the two-fermion wave function Sip^cu) 
will have two terms. The first one is the free solution (g — 0), which repre- 
sents an incident wave. The second term is the non-trivial part, arising when 
g 7^ 0, which couples different momenta and is associated with the scattered 
wave [42, 43]. Therefore, 

H^(k,q; W ) = a%-k) + (48) 

1 -32vr# 2 Ikl 2 r Ikl . „ . 

5y?(k,q;a;) , 



[u 2 - 4k 2 + ie] [u 2 - /i 2 + Zw 2 - S] A; A' A; 

where q is the relative momentum for two incident quasi-fermions and a is an 
overall phase factor. The outgoing-wave (-Me) boundary condition was used 
to solve Eq. (46), but other conditions, e. g., the incoming- wave condition 
{—ie) or Van Kampen wave condition [44], could alternatively have been 
chosen. 

The integral equation (48) can be solved as usual [43]. We integrate both 
sides with respect to k to obtain the expression 

ex 
v(k) 5<p(k, q; u) = ——a r , (49) 

u 2 -^ 2 + Zu 2 -Y,) K ' 
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where 



/+M 



kfc u 2 - 4k 2 + ie 



(50) 



with k = \/k 2 + fh 2 , while fh = m — g(0)| cq stands for the effective mass 
of fermion particle and (0)| oq is given by Eq. (31). Substitution of this last 
result in Eq. (48) yields the equality 

l k l x /, x x , n / «^0 \ |k| 1 |k| 

V 1 dipik, q;cu) = ad q- k) - — -^ - T - A _ L . , ^- (51) 

with 

A + (^) = ^— 2 (co 2 - /i 2 + Za; 2 - E) + I+(uj) . (52) 

Finally, substitution of Eq. (49) in (44) determines the oscillation fre- 
quencies 



w = 2g 2 = 2yjq 2 + fh? , (53) 

where q is the relative momentum for two incident quasi-fermions with mass 
fh. 

We have therefore found an analytical solution for the elastic channel of 
the two-fermion scattering problem defined by Eqs. (43) and (44). 

Thanks to the special form of the interacting potential, the following 
closed expression for scattering matrix can also be obtained [43] 

T(k,k';u) = v(k)^- ) v(k') (54) 

with A + (u) given by Eq. (52). 

In summary, this section discussed the solutions of the RPA equations. 
These elementary excitations describe a coupled channel scattering problem. 
The particular case of two-fermion elastic process was studied. Given the 
simple interacting potential, we were able to obtain closed expressions for 
the two-fermion wave function and the scattering matrix. Several dynamical 
behaviors can be read off from A + (w). The remaining problem is the diver- 
gent integral I + {oo) in Eq. (52), to removed with the help of counterterms. 

In Sec. 4 we will see that, in addition to the counterterms used in the 
stationary-state calculation [9], a convenient wave-function renormalization 
constant Z will have to be chosen. 
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4 Renormalization and bound state solution 



We now use the framework developed in Refs. [11, 12, 42] to investigate the 
conditions for the existence of bound states of Dirac spin-1/2 particles in a 
Higgs- Yukawa system. The standard procedure is to analyze the positions of 
the poles of the scattering matrix (54). Equation (52), however, contains a 
divergent integral. We will next show that the divergent terms can be kept 
directly under control with the help of Eqs. (21)-(24) and a convenient choice 
for Z, which yields a finite expression for A + (u). 

We therefore substitute the counterterms (21)- (24) in Eq. (52), and after 
some algebra, obtain the expression 

A+(w) = —^[(l + Z)u 2 -iJ 2 + lQ7rg 2 G(0)-247rg 2 M 2 L(m)}+I + (uj) (55) 
32ng 2 

with I + ((jj) given by Eq. (50). 

In the interval < uj < 2m the integral I u is well defined, and we can let 
e = 0. For uj > 2m, on the other hand, the spectrum defines a continuum. 
Straightfoward calculation yields the result 
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(56) 
(57) 



where a regularizing momentum cutoff A was introduced, and the finite term 
F(u) is given by the relation 

2(4m 2 - to 2 ) 3 / 2 
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(58) 
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(59) 



In Eq. (56), 9 is the Heaviside function, defined by the relations 

= if 
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From Eqs. (55)-(57), we can immediately see that there is still a logarithmic 
divergence. To cancel it, we choose the following wave-function renormaliza- 
tion [23]: 

Z = A7vg 2 L(m) . (61) 

The resulting finite expression is 



,2 



A + (uj) = -^— + F(u)-9(uj-Arri 2 )—u 2 -Am z . (62) 

g 2 m 2 8n L J 

The derivation of Eq. (62) fixed several counterterms, given by Eqs. (21)- 
(24), and (61), to eliminate the divergences. These counterterms are not 
unique, since they are well defined except for a finite value. This makes the 
results dependent on the renormalization scheme. The arbitrary finite con- 
stants are usually determined by high-energy experiments. The dependence 
on the renormalization scale is therefore often used to estimate the accu- 
racy of the theory. In the case of a scalar plasma system, or Higgs- Yukawa 
system, Refs. [9, 22, 23] have discussed the determination of these arbitrary 
finite constants to obtain the finite stationary Eqs. (26)-(31) in canonical 
form. 

We next face the problem of obtaining the poles of the scattering matrix 
when 

A + (w) = 0. (63) 

Depending on u, the system has different dynamical behaviors [42]. For 
u 2 < 0, the system is unstable, since the exponentials on the right-hand 
sides of Eq. (42) become real. For u 2 > 0, by contrast, the system is in the 
scattering regime. The solution of interest lies in the interval of < u 2 < 
4m 2 . In this interval the system may have a stable bound state if there exist 
uj b such that A + (u B ) = 0. Figure 1 shows A + (w) as a function of u/m, 
when g — 1, for three combinations of //, in unit of m. Also for g — 1, 
A + (u) has a single (no) zero when fi/m < 1.794 (fi/m > 1.794). A natural 
interpretation considers that, at fixed coupling, the boson mass determines 
the range of the Yukawa potential. When /z is large, it is more difficult for the 
fermions to interact and, consequently, the probability of forming a bound 
state, decreases. This behavior is, however, compensated by increases in g, 
as shown by Fig. 2, which plots the condition (63) in the (/x/m , g) plane. 

Figure 2 shows that the behavior of fi/m x g becomes nearly linear for 
g > 50. We can apply these results, obtained from the Higgs- Yukawa model 
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in the RPA mean-field approximation, to recent observations at the LHC 
(ATLAS and CMS collaborations), which led to the announcement of a pos- 
sible observation of a Higgs boson [13]- [15]. In the experiments five de- 
cay channels of the Higgs boson were observed, i. e., — > 77, — y bb, 
— y t + t~, — y WW and — y ZZ. The production of Higgs bosons in 
proton-proton collisions is known to occur through multiple channels, with 
branching ratios dependent on the mass of the Higgs boson. Reference [45] 
presents the branching ratios of the Higgs decay channels as a function of its 
mass. For Higgs masses below 130 GeV/c 2 , the Higgs boson is expected to 
decay mainly in the following fermions: bottom quarks b, charmed quarks c, 
and tau leptons r. 

From Fig. 2 we can determine limits to the intensity of the coupling 
constant g of the Higgs- Yukawa model, in the RPA mean-field approximation, 
for each decay channel of the Higgs boson. Let us consider a Higgs boson 
mass of m^=125 GeV/c 2 , and masses mb=4.2 GeV/c 2 , m T =1.8 GeV/c 2 , and 
m c =1.3 GeV/c 2 for the botton quarks b, tau leptons r, and charmed quarks c, 
respectively [45] . We then have the following ratios between the Higgs-boson 
mass and the masses of the three fermions: 

^-30 , ^-70 , ^95. (64) 

nib m T m c 

Therefore, for these mass ratios and no bound states of fermions (decay 
channels), the intensity of the Higgs- Yukawa coupling for the decay channel 
— y bb is limited by the condition g(<f), b) < 570. Similarly, for the decay 
channel — y t + t~ one obtain g((f>,T) < 1300, and for the decay channel 
—y cc one obtains g(0, c) < 1800. 

Figure 2 shows that, in our approach, only two-fermion bound states can 
exist for fi/m < 1.1. These results contradict the predictions of the Higgs 
boson decaying into top quarks, since m t = 173GeV/c 2 yields m^jm t = 
125/173 ~ 0.7. Our calculations indicate that only Higgs bosons with masses 
larger than 190 GeV/c 2 can decay into top quarks. These results are consis- 
tent with the results in the literature [45]. 

Finally, we want to emphasize that the phase diagram in Fig. 2 was ob- 
tained in RPA mean-field approximation. As discussed in Refs. [39, 46, 47, 
48] , the contribution of the collisional effects grows with the coupling constant 
g. For large g one finds large non-unitary contributions from the collisional 
effects. For coupling constants in the range < g < 100, the corrections, 
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i. e., the collisional terms, cannot be neglected. Systematic corrections adding 
dynamical correlation effects to the RPA mean-field calculations can in prin- 
ciple be readily obtained with the help of a projection technique discussed in 
Refs. [9, 37, 39]. The resulting occupation numbers are no longer constant 
and affect the effective dynamics of the Gaussian observables. In particular, a 
finite matter-density calculation beyond the mean-field approximation would 
allow study of such collisional observables as the transport coefficients. 

5 Conclusions 

Reference [9] treated the initial-values problem in a quantum-field theory 
of interacting fermion-scalar field theories in the Gaussian approximation. 
Although quite general, the procedure was implemented for the vacuum of 
an uniform (3+1) dimensional relativistic quantum Higgs- Yukawa model. 
The TDHF renormalized kinetic equations describing the effective dynamics 
of the Gaussian observables in the mean-field approximation were obtained. 

The present work has adapted a non-perturbative framework, the Kerman- 
Lin procedure [10-11], to investigate the near equilibrium dynamics close to 
the stationary solution of arbitrary interacting fermion-scalar field theories. 
As a application, we have chosen to describe the RPA-excitation of the Higgs- 
Yukawa system at zero temperature. 

We have studied the linearized form of the mean-field kinetic equations in 
Ref. [9] around the stationary (vacuum) solution. In this context, the RPA 
oscillation amplitudes of excitations were identified with the wave functions of 
quantum particles and the resulting equations enabled us to study scattering 
processes, non-perturbatively. These RPA equations were solved analytically 
by well-known scattering-theory procedures, which yielded a simple form for 
the scattering amplitude. We have also shown that the usual definitions of 
counterterms can be applied to the resulting expression, from which relevant 
physical aspects of the system excitations can be obtained. In particular, 
the results indicate that bound states exist in certain region of the phase 
diagram. 

We have applied our results to recent observations at the LHC ATLAS and 
CMS collaborations. We have obtained limits for the intensity of the coupling 
constant g of the Higgs- Yukawa model, in the RPA mean-field approximation, 
for three decay channels of the Higgs boson. 
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Finally, we comment that, in principle, systematic corrections can read- 
ily be added to the RPA mean-field treatment with the help of a projection 
technique discussed in Refs. [9, 37, 39]. The no-longer constant occupation 
numbers will affect the effective dynamics of the Gaussian observables. The 
framework in this paper alse serves as groundwork to discussions of finite 
densities and finite temperatures [40]. In particular, a finite matter-density 
calculation beyond the mean-field approximation allows one to study col- 
lisional observables, such as transport coefficients. The extension of this 
procedure to explore nonuniform systems is straightforward; unfortunately, 
it is tedious. In this case, the spatial dependence of the field are expanded 
in natural orbitals of the extended one-body density. These orbitals can be 
expressed in terms of a momentum expansion by means of a more general 
Bogoliubov transformation [31, 33]. 
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Figure Captions 



Figure 1. The behavior of the function A + (u) as a function of energy u, in 
unit of m, for several values of the fi/m, when g — 1 is fixed. 

Figure 2. Existence of bound state of two fermion as a function of parameters 
n/m and g. 
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